home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / DPPFA.z / DPPFA
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPPPPPFFFFAAAA((((3333FFFF))))                                                            DDDDPPPPPPPPFFFFAAAA((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPPFA   - DPPFA factors a double precision symmetric positive definite
  10.      matrix stored in packed form.
  11.  
  12.      DPPFA is usually called by DPPCO, but it can be called directly with a
  13.      saving in time if  RCOND  is not needed.  (time for DPPCO) = (1 +
  14.      18/N)*(time for DPPFA) .
  15.  
  16.  
  17. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  18.       SUBROUTINE DPPFA(AP,N,INFO)
  19.  
  20. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  21.      On Entry
  22.  
  23.      AAAAPPPP DOUBLE PRECISION (N*(N+1)/2)
  24.         the packed form of a symmetric matrix  A .  The
  25.         columns of the upper triangle are stored sequentially
  26.         in a one-dimensional array of length  N*(N+1)/2 .
  27.         See comments below for details.
  28.  
  29.      NNNN INTEGER
  30.         the order of the matrix  A .  On Return
  31.  
  32.      AAAAPPPP an upper triangular matrix  R , stored in packed
  33.         form, so that  A = TRANS(R)*R .
  34.  
  35.      IIIINNNNFFFFOOOO INTEGER
  36.         = 0  for normal return.
  37.         = K  if the leading minor of order  K  is not
  38.         positive definite.  Packed Storage The following program segment will
  39.      pack the upper triangle of a symmetric matrix.
  40.         K = 0
  41.         DO 20 J = 1, N
  42.         DO 10 I = 1, J
  43.         K = K + 1
  44.         AP(K) = A(I,J)
  45.         10    CONTINUE
  46.         20 CONTINUE LINPACK.  This version dated 08/14/78 .  Cleve Moler,
  47.      University of New Mexico, Argonne National Lab.  Subroutines and
  48.      Functions BLAS DDOT Fortran DSQRT
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.